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1. Introduction 



For at least the next decade, the energy frontier for accelerator-based particle physics will 
be located at hadron colliders, the Tevatron at Fermilab and the Large Hadron Collider 
at CERN. At a given large momentum transfer, the most copious events at these colliders 
should be hadronic jets. To test the Standard Model at the shortest possible distances, 
therefore, the jet production cross section should be known with the highest possible preci- 
sion. Existing calculations of jet production at next-to-leading order (NLO) in the strong 
coupling constant [||, ^, ^ agree well with the data over a broad range of transverse 
momentum. Still, the NLO predictions have an uncertainty from higher order corrections, 
traditionally estimated from dependence on the renormalization and factorization scales, 
which is of order 10% or more. For very large momentum transfer the predictions can be 
improved by resumming threshold logarithms Q. There are also sizable uncertainties as- 
sociated with the experimental input to the parton distribution functions [^, even though 
global fits to the data have recently been performed |^] within an approximate next-to- 
next-to-leading order (NNLO) framework 0]. Nevertheless, an exact NNLO computation 
of jet production rates would be very welcome. Besides reducing the scale uncertainties for 
jet rates, the same numerical program should allow a better understanding of energy flows 
within jets, as a jet may consist of up to three partons at this order. 

Several types of QCD amplitudes are required for a NNLO calculation of jet production 
at hadron colliders. Both the tree amplitudes for six external partons ^ and the one-loop 
amplitudes for five external partons |jl^ have been known for some time now. Recently, in 
a tour de force series of calculations, Anastasiou, Glover, Oleari, and Tejeda-Yeomans have 
provided the NNLO interferences of the two-loop amplitudes with the tree amplitudes, for 
all QCD four-parton processes, summed over all external helicities and colors ijlTI , [T^ . 

In this paper, we compute the gg gg amplitudes directly at two loops in the spinor 



helicity formalism |13], and expose their full dependence on external colors as well. The 
additional helicity and color information provided here is not necessary for the main phe- 
nomenological application, NNLO jet production in collisions of unpolarized hadrons. How- 
ever, it still provides several benefits: 

• Jet production in collisions of polarized protons, as planned for the relativistic heavy 
ion collider (RHIC) at Brookhaven, may help to determine the poorly-known po- 



larized gluon distribution in the proton |14]. Theoretical predictions of the relevant 



observables require scattering amplitudes for polarized partons. Currently, predic- 
tions are available through NLO |15|; the helicity amplitudes presented here are a 



prerequisite for improving the predictions to NNLO accuracy. 

Many formal properties of scattering amplitudes are simpler in a helicity basis and/or 
after color decomposition. Such properties include supersymmetry Ward identi- 
ties dl^], collinear limits 17, 18 1, and high-energy behavior [|l9|] . 



Our results serve as a check of the results of ref. [12|, and are useful for investigating 
the dependence of two-loop amplitudes on the variant of dimensional regularization 
used. 
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Here we also present the helicity amplitudes for gg gg scattering in pure = 1 
supersymmetric SU {N) gauge theory. Such amplitudes only differ from QCD with massless 
quarks in that the fermions are in the adjoint rather than the fundamental representation; 
yet they obey supersymmetry Ward identities p^ ] and are generally simpler than their 
QCD counterparts. They also provide useful auxiliary functions for describing the QCD 
results. 

Several versions of dimensional regularization have been used for loop calculations in 
QCD, differing mainly in the number of gluon polarization states they assign in 4 — 2e 
dimensions. The conventional dimensional regularization (CDR) scheme [20| assigns D — 
2 = 2 — 2e states to all gluons, whether internal or external, virtual or real. This scheme is 
traditionally employed in calculations of amplitude interferences, such as ref. []T2| . In the 
helicity approach, the number of external, observed gluon states is necessarily 2 (helicity 
±1), but there is some freedom in the number of virtual gluon polarizations. The 't Hooft- 
Veltman (HV) scheme |21] contains 2 — 2e virtual gluon states, while the four-dimensional 
hehcity (FDH) scheme [||, |3| assigns 2. The FDH scheme is related to dimensional 
reduction (DR) |^4| but is more compatible with the helicity method, because it allows 2 
transverse dimensions in which to define helicity. Of these variants, only the FDH scheme 
is fully compatible with supersymmetry Ward identities for helicity amplitudes, some of 
which have been verified through two loops ||2^. Here we work primarily in the 't Hooft- 
Veltman (HV) variant of dimensional regularization , but we also discuss the conversion 
to the CDR and FDH schemes. 

Two-loop scattering amplitudes in massless QCD possess strong infrared (soft and 
collinear) divergences. Using dimensional regularization with D = 4 — 2e, the amplitudes 
generically contain poles in e up to However, these divergences have been organized 

by Catani |25| into a relatively simple form, which is completely predictable through at 
least order l/e^. We shall use Catani's formulae and color space notation to organize the 
gg gg helicity amplitudes into singular terms (which do contain terms in their series 
expansion in e), plus finite remainders. We find that the general form of the divergences 
given in ref. |2^ holds precisely in both the HV and FDH schemes; however, the numerical 
value of the coefficient K, which appears at order differs in the FDH scheme from its 
value 1 25] in the HV (or MS) scheme. 

The 1/e poles were not predicted a priori in ref. [^] for general processes at two loops. 
For the gg gg amplitude, ref. computed the interference of the 1/e pole terms with 
the tree amplitude, summed over all colors and helicities. Here we extract the full color 
and helicity dependence of the 1/e pole terms. We find a term which is independent of 
color and helicity, and which agrees with that found by ref. |jl2| (when we use the HV 
scheme), plus a second term with nontrivial color-dependence, which vanishes when the 
color-summed interference is performed. A term with similar color structure has also been 
identified in contributions of one-loop factors for soft radiation to NNLO processes [ p6| . 
We shall also discuss how terms in the infrared decomposition of ref. ||2^ are modified, 
beginning at order 1/e^, in other variants of dimensional regularization, such as the FDH 
scheme. 

The paper is organized as follows. In section we review the infrared and color struc- 
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ture of one- and two-loop QCD amplitudes. In section ^ we describe the one-loop gg —>■ gg 
amplitudes in a form that is valid to all orders in e |17, 27, 28, and show how to expand 
them through 0{e^). This accuracy is required because one- loop amplitudes enter the for- 
mulae for the singular parts of two- loop amplitudes multiplied by . Section shows 
that apart from this requirement, only finite remainder terms in the one- loop amplitudes 
are needed, because of cancellations with other NNLO contributions. These remainder 



terms are then tabulated in section 3.3 



In section]^ we describe our method for computing the two-loop amplitudes. Section 4.1 



summarizes how we evaluate loop integrals, especially those that arise only in the helicity 



method. Some consistency checks on the results are listed in section [4.2| . Section |4.3| dis- 
cusses the additional singular term appearing at order 1/e in the color-decomposed gg — > gg 
amplitude, which does not contribute to the color-summed interference with the tree am- 
plitude. The finite two-loop remainder functions in the HV scheme are then presented in 
section 4.4 and appendix 

In section ^ we describe conversion of the HV results to different schemes, and the com- 
parison with ref. |T^, after our results are summed over all external colors and helicities. In 
section |6| and appendix ^ we give the two- loop amplitudes for pure = 1 supersymmetric 
Yang-Mills theory, whose finite remainders also serve as auxiliary functions for describing 
the QCD results. In section we present our conclusions. 



2. Review of infrared and color structure 

In this section we review the structure of the infrared singularities of dimensionally regu- 
larized one- and two-loop QCD amplitudes, using Catani's color space notation [^, as a 
prelude to presenting the finite remainders of the one- and two-loop gg — > gg amplitudes. 
The process considered in this paper is 

g{-pi, -Ai) + g{-p2, -A2) g{p3, A3) + 5(^4, A4) , (2.1) 

using an "all-outgoing" convention for the external momentum (pi) and helicity (A^) label- 
ing. The Mandelstam variables are s = {pi +P2)'^, t = {pi +^4)^, and u = {pi +^3)^. 

We work with ultraviolet renormalized amplitudes, and employ the MS running cou- 
pling for QCD, Os^fi). The relation between the bare coupling a" and renormalized coupling 
as(//), through two-loop order, is 

< nl' S, = a,(/z) fi'^' 

where n is the renormalization scale, = exp[e(ln47r + V'(l))], and 7 = —ip{l) = 0.5772 . . . 
is Euler's constant. The first two coefficients appearing in the beta function for QCD, or 
more generally SU{N) gauge theory with flavors of massless fundamental representation 
quarks, are 

^ UCA-4:TRNf ^ 17Cl-{WCA + 6CF)TRNf 

"0 = 7 ) "1 = 7 ' y^-^) 



2-K e 



2lT 



(2.2) 
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where Ca = N,Cf = (iV^ - l)/(2iV), and Tr = 1/2. (Note that ref. [||] uses the notation 
/?o = W(2vr),/3i = 6i/(2^)2.) 

The perturbative expansion of the gg — > gg amphtude is 



M 



99~*aa 



2 



(2.4) 



27r 



+ ( ^) -M£L,,(M;M) + 0(a?(M)) 



where ■M.'ggUggif^] {p}) is the L}'^ loop contribution. Equation (2^) is equivalent to the 
following MS renormalization prescriptions at one and two loops, 



•' 99-'9g 

•' gg-'99 



^-l^(l)unren _ ^_y^(0) 
e 99-^99 £ 99^99 ' 



(2.5) 



2e 



The infrared divergences of renormalized one- and two-loop n-point amplitudes are 
given by [25|, 



|A^«(/x; M)),s = /i; M) M))Ka + \M^^^'''' {l^; M))Ra , (2.7) 

|>l(2)(/x;M)),s = /(^)(e,/i;M) |Al«(/i;M))aa (2.8) 
+ I^ll (e, ^; M) W (/i; M)),s + l^^i'^'^'^l^; M)>aa , 

where the "ket" notation \M.\^\^\ {p})) ^s. indicates that the L-loop amplitude is treated 
as a vector in color space. The actual amplitude is extracted via 

M„(l"S...,n''") = {ai,...,an\Mn{pu...,Vn)) , 



(2.9) 



where the are color indices. The subscript aa indicates that a quantity depends on 
the choice of renormalization scheme. The divergences of Mn"^ are encoded in the color 
operator while those of Mn^ also involve the scheme-dependent operator /aa- 
In QCD, the operator I^^^ is given by 



1 p-^V-Ci) " " 



£2 + T? e 



2pi ■ Pj 



(2.10) 



where Xij = +1 if i and j are both incoming or outgoing partons, and Xij = otherwise. 
The color charge Tj = {T^} is a vector with respect to the generator label a, and an 
SU{N) matrix with respect to the color indices of the outgoing parton i. For external 
gluons = so Tf = Ca = N, and 



75 



IICa - 4.TRNf 



(2.11) 
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The operator Ij^^l is given by ||2f 

where the coefficient K^s. in either the HV or CDR schemes is given by |2 



67 



vr 



10. 



(2.12) 



(2.13) 



Although no scheme dependence was assigned to this coefficient in ref. we shall find 
in section |5| that it is scheme dependent. The function H^l contains only single poles, 



H'il{e,fi;{p}) = 0{l/e), 



(2.14) 



but is not predicted a priori for general processes. The color- and helicity-summed ma- 
trix element {M^^^\H^'^\e)\A4^^'^) has previously been computed in the CDR scheme for 
gg gg [p!2| (and for some other multi-parton processes |11, |30|]). We shall extract the full 

^ — (2) — — 

color and helicity dependence of -ffjj^g. (^) 99 ~^ 99 the HV scheme in section 4.3, and 
in the FDH scheme in section ^. 

An explicit color basis for the gg gg amplitudes is given by 



-^AiA2A3A4 ~ •S'AiA2A3A4 X X/'^^'*'-^AiA2A3A4 ' 



where 



(2.15) 



tvp: 

TrW 

iv[5: 



i=l 

ti(T'^^ T""^ T'^^ T""^ ) 
tr (T'^^ T"'^ T'^^ T"'^ ) 
tr (T""^ T"^'* T'^'^ T"''^ ) 

tr (r"^^ 2^13 2^(12 rpai ^ 

tr {T""^ T""-^ T"'^ T""^ ) 
tr {T""^ T'^'^ T"^"^ T""^ ) 
tr(T"ir"2)tj.(r''3 7^a4) ^ 

tr(r'*ir"3)tr(r'*2r"4), 
tr(T"ir"*)tr(r"2r''3) . 



(2.16) 



Here T'^ are SU{N) generators in the fundamental representation, normalized according 
to the convention typically used in helicity amplitude calculations, tr(T"T'') = 6°'^. (The 
T"" used in this color decomposition should not be confused with the appearing in 
which are in the adjoint representation; nor should they be confused with the generators 
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for the quark representation, which have the more "standard" normahzation, Tr = 1/2, as 
mentioned above.) 

We have also taken the opportunity in eq. ( |2.15| ) to remove some hehcity-dependent 
overall phases, which arise because we evaluate the amplitudes in the spinor helicity for- 
malism [in 



s. 



. [12] [34] 
'(12) (34)' 
. (12) [34] 
'[12] (34)' 



, (12) (14) [2 4] 

(3 4) (2 3) (2 4) ' 
. (13) [2 4] 
'[13] (24)- 



(2.17) 



The spinor inner products [jl^, ^ are {ij) = (i~|i^) and [ij] = where are 

massless Weyl spinors of momentum ki, labeled with the sign of the helicity. They are 
anti-symmetric, with norm | (i j) | = | [ij] \ = s/Sij^ where sij = 2ki ■ kj. It follows that the 
Sx-^^x^X3\4 are indeed phases. They will cancel out from (and therefore may be freely omit- 
ted from) all transition probabilities involving unpolarized gluons, or circularly polarized 
gluons. 

In the basis ( |2.16| ) for gg gg, the matrix I^^^ is [^] 



r(i 



1 _^ 

e2 + Ne 



(2.18) 



/iV(SH 


-T) 



























(T - U) 


(S-U)\ 







iV(SH 


HU) 






















(U-T) (S-T) 















N{T + U) 



















(T-S) 


(U-S) 















7V(TH 


-u) 














(T-S) 


(U-S) 




















7V(S + U) 









(U-T) (S-T) 






























hT) 


(T - U) 


(S-U) 


(s- 


u) 


(s- 


T) 










(s 


-T) 


(s 




u) 


2iVS 










(u- 


T) 


(U-S) 


(u- 


s) 


(u 


-T) 









2iVU 





V (T- 


u) 







(T-S) 


(T- 


s) 







(T 




u) 





2m ) 



where 



2 \ <^ 



u 



(2.19) 



A reflection identity implies that the coefficients of two color structures with reversed 
j'o.i ordering are identical, so that 



M 



A1A2A3A4 



M 



(i)-[3] 
Ai A2 A3A4 ' 



^(i),[5] 
A1A2A3A4 



^(i),[2] 
Ai A2A3A4 ' 



M 



{L)\Q\ 
Ai A2A3A4 



M 



Ai A2A3A4 ■ 



(2.20) 



Also, due to Bose symmetry, parity, and time-reversal symmetry for the process ( |2.1| ), we 
only have to give results for the four helicity configurations 



A1A2A3A4 = ++++, -+++, 
The tree amplitudes are given in the basis ( ^.16 ) by 



-++, 



-+-+. 



(2.21) 



M 



(o),H 



0, 



for all z. 
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M 



(0),[7] 
Ai A2A3A4 

(0),[1] 
++ 

(0),[1] 



M 



(0),[8] 



M 
M 



Ai A2A3A4 
S 
t 

u 

St 



M 



^(0),[9] __ 
Ai A2 A3A4 

(0),[2] _ S 



0, 



M 



{0),[3] 



for all Aj, 

,2 



— , M_ 



(0),[2] 



s 



^^(0),[3] 



s 
tu ' 

M 



(2.22) 



A typical partonic cross section requires an amplitude interference, summed over all 



external colors. Such interferences are evaluated in the color basis (2.16) as 

9 



(L,L') ^,j^(L) \^{L') 
A1A2A3A4 \ A1A2A3A4I A1A2A3A4'' 



where the symmetric matrix CCij = ^ 



colors 



^ ^ "^Ai A2A3A4^^*i'^Ai A2A3A4 ' 
TrI'l *Tr[-''] is 



(2.23) 



12 



/ Cl C2 


C2 


C2 


C2 


C3 


NV 


-N 


NV \ 


C2 Cl 


C2 


C2 


C3 


C2 


NV 


NV 


-N 


C2 C2 


Cl 


C3 


C2 


C2 


-N 


NV 


NV 


C2 C2 


C3 


Cl 


C2 


C2 


-N 


NV 


NV 


C2 C3 


C2 


C2 


Cl 


C2 


NV 


NV 


-N 


C3 C2 


C2 


C2 


C2 


Cl 


NV 


-N 


NV 


NV NV 


-N 


-N NV NV iVV 


Ar2 


Ar2 


-N NV NV NV NV 


-N 


Ar2 


iVV 


Ar2 


\NV -N NV NV 


-N NV 


Ar2 


Ar2 


N^V ) 



(2.24) 



with 

Cl = iV^ - 3A^^ + 3, C2 = 3-iV2^ C3 = 3 + iV2^ V = N"^ 
The unpolarized partonic cross section is obtained from the helicity sum 

J(L,L') ^ V- j(L, 



1. 



r(L,L') 
1 A2A3A4 



(2.25) 



(2.26) 



Ai=±l 



after the usual averaging over initial spins and inclusion of flux factors. For example, the 
helicity sum for the tree-level cross section, constructed from eq. ( |2.22| ) in either the HV 
or FDH scheme, is 



j(o,o) ^ igjvVl 3 



tu 

12 



US 



St 

-7/2 



(2.27) 



3. One-loop amplitudes 



The one- loop amplitudes for gg — > gg were first evaluated through 0{e^) as an interference 



with the tree amplitude in the CDR scheme |32]. Later they were evaluated as helicity 



amplitudes in the HV and FDH schemes |3^ . 

Because /^"^^ contains terms of order 1/e^, the I^^^\Ai^^'^)^,s. term in the infrared de- 
composition (2^) of the two- loop gg gg amplitudes requires the series expansions of the 
one-loop amplitudes through ©(e^). In ref. [^], using results from refs. |17, 27, 25], the 
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one-loop gg gg helicity amplitudes were presented in a representation valid to all orders 
in e, in both the HV and FDH schemes. These results can easily be rewritten in terms 
of integral functions whose series expansions are known to the requisite order [34, In 



section 3.1 we present the all-order results in the color basis ( |2.16| ), with the normalizations 



implicit in eq. (p.4|). 

In section we show that the only place that terms beyond 0{e^) in the one-loop 
amplitudes are required in an NNLO calculation is in the infrared decomposition ( |2.8D of 
the two-loop amplitudes. 

Finally, in section |3.3| we list the finite remainders of the one-loop amplitudes in the HV 
scheme, after the renormalization (p.5D and subtraction of infrared divergences (2.7). The 



corresponding finite remainder in the one- loop/one- loop NNLO interference has already 
been computed in the CDR scheme, summed over all colors and helicities [^]. Our HV 
amplitude remainders lead to precisely the same result. 

3.1 All orders in e 



Here we present the renormalized one- loop gg gg amplitudes in the color basis ( 2.1(^ ), 
with the normalizations implicit in eq. (p.4|), in a form valid to all orders in e. 



The first coefficient in the color basis ( |2.16 ) for gg gg at one loop may be written 



in terms of "primitive" amplitudes for a gluon or quark in the loop, as |36| 

Mi'^lxMt,u) = + NfMi—^ - 7<S3A. ■ (3.1) 

The remaining single-trace coefficients are obtained via crossing symmetry: 





t, 


u) 


— m^^^j^ 


{s, 


u,t), 




(s, t, u) 


- M^'^'W 


(n. 


t,s) 




t, 


u) 


= M « ' W 


(s, 


u,t), 




(s, t, u) 


- M^'^'W 


(u, 


t,s) 




t, 


u) 


= M«' W 


(s, 


u,t), 




(s, t, u) 


= M« ■ W 


(m, 


t,s) 




t, 


u) 


= M«' W 


{u 


s,t), 




{s,t,u) 


= M« 


(n, 


t,s) 



(3.2) 

where appropriate analytic continuations are required to bring each function into the phys- 
ical region. The double trace coefficients, to which only the gluon loops contribute, follow 



from a [/(I) decoupling identity |36]: 



j^mr] ^^(1),[8] ^^^(1),[9] ,Af(')J2] ,^j{l).[3] 

A1A2A3A4 A1A2A3A4 A1A2A3A4 jY V A1A2A3A4 A1A2A3A4 A1A2A3A4 y ^ _Q 

luon '^^'^^ 

It is convenient to write the gluon and fermion loop contributions, Aff|A°A3A4 
-^Ai'a2A3A4' ™ terms of a supersymmetric decomposition into scalar, chiral = 1, and 
= 4 supersymmetric multiplets in the loop |l^: 

<rA3A4 = (1 - ^^n)MZtl,x. - ^Kxkx. + <r'A3A4 > (3-4) 

J ^-fermion _ _ ii ,rscalar l\/fN=l fi ^\ 

^^^^1X2X3X4 - ^'^AiA2A3A4 + ^'^AiA2A3A4 ' l-^-^J 
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where = 1 for the HV scheme and Sr = for the FDH scheme. 

For "maximally helicity violating" configurations, the supersymmetric components 
vanish by a supersymmetry Ward identity fl^ , 



M^Vt+ = M^^l^ = M^=t^ = M^=t^ = 0. (3.6) 



The remaining independent components are [17, 27, 29| 

^scalar^ = -e(l - e) BoX^^\s, t) , 

^scalar ^ ^(^i^, Xri(6) (.) + fi^i^e Tri^^) (t) + 
su tu 

+ izJ^eBub(6)(s) + ^eBub(6)(t) - ^eBox(6)(s,t) - 

2ii 

-e(l-e)Box(^)(s,t), 

^scalar^ = (t) - e{l - e) Box(8)(s,t) , 

™iar ^ _i^g^i(6)(^) _ i^gTri(6)(s) + -Bub(4)(t) + -Bub^^^fs) + 
u u u u 

+ !i^Bub(^)(t) + ^i^Bub(6)(s) - Tri(6)(t) - Tri(6)(s) - 

- -Box(6) (s, t) - e(l - e) Box^^) (s, i) , (3.7) 
u 





1 

~ ~2 




~ 2 s 




1 . 


m5-^ 


1 

= 2" 



' 2t 

i^But 
2 t 



(3.8) 



)x(^)(s,t). (3.9) 

Here Bub('')(s), Tri('^)(s) and Box(")(s,t) are the one-loop bubble, triangle and box scalar 
integrals, evaluated in D = n — 2e dimensions. The bubble and box integrals are 

= - 2.(1 -20(3 -2.) 



where 



r(l + e)r2(l-e) 
r(l-2e) 



l-lc2e'-lCse'-%ae' + 0{e'), (3.11) 
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with 



n 



C2 



n=l 



vr 
~6 



C3 = 1.202057... 



Ca 



TT 

90 



(3.12) 



and we have kept the full dependence on e in the integrals. In the s-channel (s > 0), 
e-expansions of the functions ( |3.10| ) are given by using the analytic continuation ln(— s) — > 

In S — ITT. 

The box integrals in various dimensions appearing in eqs. (|3.7| ) and (|3.9| ) are related 
via a dimension-shifting formula valid to all orders in e, 

1 



Box(s)(s,t) 



stBox(4)(s,t) - 2tTti^^\t) - 2sTri(4)(s 
St Box(^) (s, t) - 2t Tri^^) (t) - 2s Tri^^) (s) ) . 



(3.13) 



2(-l + 2e)n 
1 

2(-3 + 2e)n 

Because the D = 6 — 2e scalar box integral is completely finite as e — > 0, it is convenient 
to express the other box integrals in terms of it. This isolates all divergences to triangle 
and bubble integrals. To expand the six-dimensional box to higher orders in e, one could 
use an expression for Box*-'^^(s, t) valid to all orders in e, in terms of hypergeometric func- 
tions [39, |3^, and the dimension-shifting formula ( 3.13| ). Or one can expand the Feynman 
parameter integrand for Box^^\s,t) in e directly. 

In the u-channel (s < 0, t < 0), where the functions are manifestly real, the expansion 
of the six-dimensional box through 0{e^) is [34, 35] 



Box(6)(s,t) 



rru 



-l-e 



2(1 - 2e) 



( y3 ^2 ^ 

+ 2e hh{-v) - VU2{-v) - ^ - 



2e2 ( hu{-v) + WU^{-v) - -V'^U2{-v) 



7r2 



1 



1 



vr 



6 



- — VW-2CA\+{s^t) 



where 



w 



u 



4 



V = \n\^--j, W = ln[ 
In the s-channel (s > 0, t < 0) an analytic continuation of the box integral yields, 
Box(6)(.,i) "r[^[-^ ,, 



(3.14) 



(3.15) 



m(1 - 2e) 



U3{-x) + XU2 
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1 



Li4(-y) +Li3(-y)X + 



+ ivr 



+ ix3y - i(x2 + vr2)2 + ^XY + —A + 
3 8^ ^3 360 y 

X + e(Li2(-x)+yX-ix2-^) + 

+ (-Li3(-x) - Li3(-y) - \yx'' + ^-X^ + Cs 



+ 



+ 0(e3) 



where 



X = In 



y = In 



Also define, for future use, 



X = X + m, 



Y = Y + i7r. 



(3.16) 
(3.17) 

(3.18) 



In general, both expansions of box integrals ( |3.14| ) and (3.16) will appear in eq. ( p.8| ) for 
the divergences of the two-loop amplitudes. 

3.2 NNLO cancellations involving M.^gg^gg 

The NNLO gg gg cross section has a term proportional to the square of the one-loop 
amplitude, /(^'i) = {M.^g^^„g\J^^ilfj_^gg) . One might expect to need the 0{e'^) terms in 
M.^gg^gg here, because M.gg-*gg also contains l/e^ terms. Here we show this is not the 
case, for a generic NNLO calculation of an infrared-safe observable, because of cancellations 
against contributions involving radiation of additional partons. 

If one uses the one- loop infrared decomposition ( |2.7| ), one can rewrite /(I'l) as [P 



J(i,i) = (Al_W_^,J/«t(,)j(i)(,)|^^£)_^^^) + 2Re(H^,)fi^^^ 



99^99 

/7U(l)fin I M^^)^'' ) 
\ 99^99' 99^99' 



99^99' 



99^99' 



+ 



-(A^SL,J/«n^)/(^n^)|A^SU,) +2Re(A^£L,J/W(e^ + 
-{M 



(l)fin |A//(l)fiii \ 
99^99'-^ 99^99' 



(3.19) 



Similarly, the contribution of the two-loop/tree interference to the NNLO gg gg cross 
section is 



2Rel(2'0) = 2Re{M'~S^gg\I^'\e)\M^^^^gg) + 2Re{M'^^^^gg\^^^^^ + 



99^99'-' 99^99' " 

(For clarity, we have dropped as subscripts from equations in this subsection.) 



(3.20) 
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Thus the two singular terms containing the one-loop amplitude in the NNLO gg gg 
cross section combine to form 

2Re(A^g_.,J(l(i)(6) + /(^)t(e))|A^£L,,) ■ (3.21) 

These terms will partially cancel, in an NNLO cross section for an infrared-safe quantity, 
against phase-space integration of certain terms arising from the one-loop / tree interference 
for the processes with one additional parton radiated (in this case, gg — > ggg, gq ggq, 
gg — > gqq, etc.). The "radiation" terms may be written generically as 

2Re{M^}jMflj . (3.22) 



To see this NNLO cancellation, it is useful to recall the corresponding cancellation at 
NLO, where the singular part of the virtual correction 



Re{M^S^J{I^'\e) + /(^)t(,))|^(o)_^^^) (3.23) 
is cancelled by phase-space integration of the real radiation terms, 

Re(A^£Ul-^Sad) ■ (3-24) 

The singular phase-space behavior, soft or collinear, of the one-loop five-point amplitude 
factorizes as 0, H, ||, |l| 



^i]L -^I'^'^W + M^S^'^ , (3.25) 

where {S'-^^) represents a universal tree-level (one-loop) soft or collinear factor, which 
contains all the dependence on the unresolved phase-space variables that have to be inte- 
grated over. The tree-level factorization is of course 

Mfl, Mf^S('^ . (3.26) 



Using eqs. ( 3.25| ) and ( 3.26| ), one sees that the terms in the singular behavior 

of eq. ( p. 22 ) have exactly the same form as the real NLO terms ( |3.24| ), but with {M^^^^ 
<(i)| 



replaced by \ , and an overall factor of 2 from the interference. Thus the result of 

integrating the M^^S^^^ terms in eq. ( 3.22| ) over phase space must be cancelled by the 



virtual NLO terms ( p.23| ), but with the corresponding replacements, i.e. by eq. (3.21). 



The NLO cancellation is good to 0{e^) (after factorizing initial-state collinear singularities 
in the usual manner); see e.g. ref. [^2[. The NNLO cancellation is at the same order, in 
the sense that 0{e^) and higher terms in Ai^^^ no longer contribute. 

In summary, the only place the terms beyond ©(e*^) in A4^^ are really required at 
NNLO is in the infrared decomposition of the two-loop amplitude. Once the two-loop 
finite remainders A4^^^^^ are given, the higher-order terms in Ml^^^ are no longer needed. 
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3.3 Finite remainders 

Next we tabulate the finite remainders of the one-loop gg —>■ gg amplitudes at 0{e^), defined 

^(l),[i]fin 
A1A2A3A4 



by M-^^g^gg 111 eq. ( |2.7| ) and color decomposed into -^{J^;^2*A^^4 in eq. ( 2.15| ). We write, in 



the HV scheme, 

<S:x, = -bo{Hs/^^') - ^^)MrZ.^. + ^ «Sa.A3A. + ,,,3,, , ^ = 1,2,3, 

-^Aia1a3A4 = 9a1a2A3A4 + ^a1a2A3A4 ' ^ = 7, 8, 9, (3.27) 

and the M^^^'W'^'^ for i = 4,5,6 follow from eq. ( p. 20 ). The one-loop ^7(1) decoupling 
identity ( |3.3P implies that 

5'a1a2A3A4 ^ ^('^AiA2A3A4 + "'AiA2A3A4 + °'AiA2A3A4) ' i = 7, 8, 9. (3.28) 

For the -|— 1-, — and 1-+ helicity configurations, Bose symmetry under exchange 

of legs 3 and 4 (t n) implies that 

41a2A3A4(^'^''") = «a]a2A3A4(«'^'*)' 

41a2A3A4(^'*'") = 4a2A3A4(«'^'*)- (3-29) 

For the -|--|— helicity amplitude, the independent remainder functions a, c, g and h 



are 



For — they are 



"++ 


1 




f+ - -g> 




[3] 


1 




f+ - -g' 






1 

f+ - g. 




^++ 


1 

f+ - g> 




''■++- 


f+ ~ "■++++ 


- /it^i - 


"-+- 










[3] 


1 

'^'^ Qxy ' 














^-+- 


1 




Qxy ' 




"--+- 




- /I'^l - 



(3.30) 

(3.31) 

(3.32) 

(3.33) 
(3.34) 

(3.35) 
(3.36) 

(3.37) 
(3.38) 
(3.39) 
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IS 



For they are 

aLi,,^-^(X^ + .^) + 3^, (3.40) 

--++ 4xy \^ ' ^ \Qy 2 J 36xy 4 



+ 



|t ^ u| , 



[3] a;^ / ,,,0 o\ 1 _ \ ~, 29 5 x 



lt ^ lil , 



For — I h, the required functions are 

fil (l-a;y)^v>2 2^ X IN- 67 y2 i 



Ji] _ f JL - 1.^ (^jf^ + vr^) + ^(^~^) + 5 X - (3-2y)(3x-2y) 



2y^ 4y ' ' 6y 18x 

y 

3 



"+"+ 3x 9 X ' 



3.41) 
3.42) 

3.43) 
3.44) 



3.45) 

3.46) 
3.47) 
3.48) 

3.49) 
3.50) 
3.51) 



The contribution of the one-loop finite remainders to the NNLO gg — > gg cross section 



A,=±l 



Using the color sum matrix OOij in eq. ( 2.24| ), the color and helicity sum in /(i'^)'^'^ may be 
evaluated in terms of the above explicit expressions (|3.27l) -( 3.51 ) for ^Xj^x^x-iXi' repro- 
duces precisely the finite remainder function Finite{s, t, u) given in eq. (3.22) of ref. |31| for 
the corresponding quantity evaluated in the CDR scheme; the HV/CDR scheme difference 
for /(^'i) has been completely absorbed into the first two of the three terms in eq. (3.19). 
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1^ ^4 



4 1 ^ ^4 



Figure 1: Some of the two-loop diagrams for gg gg. 



4. Two-loop amplitudes and finite remainders 

A generic sample of two-loop Feynman diagrams for gg — > gg is shown in figure |l]. How- 
ever, we did not evaluate the diagrams directly. Instead we computed the unitarity cuts 
in various channels, working to all orders in the dimensional regularization parameter 
e = (4 — L')/2 17, 15]. Essentially we followed the approach first employed at 
two loops for the pure gluon four-point amplitude with all helicities identical [^] and for 



= 4 supersymmetric amplitudes |44]. These amplitudes were simple enough that a com- 
pact expression for the integrand could be given. The fermion loop contributions with all 
helicities identical are about as simple |^]. For other helicity configurations, the integrands 
become rather complicated. We therefore used the general integral reduction algorithms 
developed for the all-massless four-point topologies in order to reduce the 

loop integrals to a minimal basis of master integrals. To efficiently incorporate polarization 
vectors of gluons with definite helicity requires some minor extensions of these techniques, 
which we now discuss. 

4.1 Tensor loop integrals 

Here we discuss techniques for evaluating the loop integrals required for the two-loop 
amplitudes for gg gg and related processes, with an emphasis on the additional types 
of integrands encountered in the helicity amplitude method. 

In calculating a typical two-loop scattering amplitude in QCD, a large number of 
two-loop integrals are encountered. The most complicated topologies are the planar and 
non-planar double box integrals, displayed in figure |2|, which are given by 

f d^p d^q V 



(27r)^ (27r)« p2 ^2 + ^y^p _ ^^)2 {p - ki - fcs)^ (g - k^f {q-k^- k^f 



(4.1) 



^-1 



(2^ (2^ p2 g2 + _ fc^)2 _ ^2)2 (p + g + k^Y {p + q + k3 + k^Y " 

(4.2) 



Here p and g are the loop momenta, and ki, i = 1,2,3,4, are the external (outgoing) 
momenta. P is a polynomial (or tensor) in the loop momenta p and q, which accompanies 
the scalar propagator factors shown in the figure. It is generated by the numerator algebra 
of the Feynman diagram, or unitarity cut, that is being evaluated. 

/3 P + q /3 





Figure 2: The planar and non-planar double box integrals. 

In the interference method, as recently applied to two-loop QED and QCD scat- 
tering amplitudes p5|, O, O, pOl, one sums over all external polarization states in D 



dimensions. In this case, V can only depend on the loop momenta, p and q, and external 
momenta, ki. By Lorentz invariance, this dependence is only through scalar products, 



^interf . 



V{p'^,p- q,q'^,p- ki,q - ki 



(4.3) 



In contrast, in the helicity amplitude method |13] used in the present paper — and 
previously applied to two-loop amplitudes in = 4 super- Yang- Mills theory ||4^, gg gg 
for identical helicities in pure Yang- Mills theory ||3^, and the QCD processes 55 — > 77 [49| 
and 77 — > 77 |^0| — V also depends on the polarization vectors for the external 
gluons. We take D > 4 in the calculation, i.e., e < 0, in order to have two transverse 
dimensions in which to define helicities. Because the polarization vectors are intrinsically 
four-dimensional, their Lorentz products with the loop momenta distinguish between the 
four-dimensional and (— 2e)-dimensional components of p and q. Write 



?[4] + ^^9' 



(4.4) 



where p[4], (^[4] are the four-dimensional components and Xp, Xq are the (— 2e)-dimensional 
components. We use the Minkowski metric with signature (1,-1,— 1,-1,...), and write 



2 _ 2 x2 

P — P[4] ~ '^p' 



f4] - Ag, {p + qf = (P[4] + q[4]? - >^l+q, 



(4.5) 



where Xp = Xp ■ Xp > 0, Xp_^_q = {Xp + Ag)^ = Xp + Xq + 2Xp ■ Xq. Then the generic 
polynomial encountered in the helicity amplitude method has the form 



Vhei. = Vip ,p-q,q ,p- ki,q- kf, X^, X^, Xp^g,ei ■p,ei ■ q). 



(4.6) 



-16- 



We rely on reduction algorithms developed to handle general tensor integrals for the 
all-massless planar [45| and non-planar |^] double boxes, and related topologies such as 
the pentabox 1^^. These algorithms were derived using integration by parts and (for 
the non-planar double box) Lorentz invariance |48| identities, which act in the space of 
integrals with V = 1, but with the scalar propagators raised to arbitrary integer powers 
z/j. For example, for the planar double box topology one considers 



where 



if (1/1,1/2,1/3, •••,1^7) = (4vr)^y" 



(2vr)« (2vr)^ l{ (pjr 



(4.7) 



Pi = q, P2 = q-h-h, P3=P, P4=p-ki-k2, 

P5 = p-h, P6=p + q, P7 = Q-h, 



(4.8) 



and Vi G {0,1,2,...}. The reduction algorithms reduce any such integral to a linear 
combination of simpler "boundary" integrals, where at least one of the i/j vanishes, plus 
one or two master integrals with the same topology. All told, there are 10 different master 
integrals for the massless 2 — > 2 processes |45, |46| ]. 

Given an integral with V of the form ( |4.3| ) or ( |4.6| ) , it is simple to convert it to integrals 
of the form (|4.7]) using Schwinger parametrization [E2|. For the scalar integrals, using 



Pi 



dtiexpi-tipj 







and performing the p and q integrals, leads to 

7 



/f(i,i,...,i) =n 



i=l 



dti 







A(r) 



-D/2 



exp 



Qxis,t,ti 
A{T) 



A{T) = TpTg + TpTpg + TgTj 



q^pq , 



(4.9) 



(4.10) 
(4.11) 



where X labels the topology of the integral (P, NP, etc.). In eq. ( 4.11 ), Tp, Tq, Tpq are 
the sums of Schwinger parameters along the lines carrying loop momenta p, q, p + q, 
respectively. Equivalently, they are Schwinger parameters for the two-loop vacuum graph 
obtained by omitting all the external lines, as shown in figure ^. For the planar double box 
integral, with propagators numbered by eq. (|4.8|), the Ta are given by 



ts + ti + t^, 



tl+t2+ t7, 



'-pq 



tQ. 



(4.12) 



The quantity Qx is more cumbersome. For the planar double box integral, its expres- 
sion is 



Qp{s,t,ti 



tlt2it3+tA + t5) + t3ti{ti + t2 + t7)+tQ{ti + t3){t2+ti) - t t^tetj . (4.13) 



However, the precise form of Qx will not matter in the following. 
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Figure 3: The two-loop vacuum graph obtained by omitting external momenta, and its three 
Schwinger parameters, which are relevant for general two- loop integrals containing v ■ p, v ■ q, Xp, 
and Xq. 



If a polynomial V is present, Schwinger parametrization converts it to a polynomial 
in the Schwinger parameters, along with inverse powers of A. Then the more general 
parametrization, 

1 1 

^^ = YT^ dtif^'-'eM-tiPl), (4.14) 

l^ij Jo 



can be used to rewrite Schwinger parameter monomials as integrals of the form ([4.7]), 
typically in shifted dimensions, D ^ D + 2n, n = 0, 1, 2, . . . (to account for the inverse 
powers of A). (Shifted-dimension integrals pose no problem; equations for them can be 
found by rewriting the factor A~^/^ in eq. ( [4.10D as 

and reducing the latter, shifted-dimension representation.) In principle, this approach gives 
a prescription to handle any polynomial in the loop momentum, for either the interference 
or helicity method. 

However, as the degree ^ • Vi increases, the number of integrals of the form (^]^) grows 
rapidly, and the reduction algorithm can become rather time-consuming. We have found 
it useful to instead use simple algebraic relations, e.g. for the planar double box integral, 

2q-k4^= p\- p^, 2q-k3=pj-pl + 2k3 ■ k^, 
2p-ki = pI- pI, 2p-k2=pl-pl + 2ki ■ k2, 
2p-q=pl-pj-pl (4.16) 

to quickly reduce integrals with polynomials of the form Vmteri. to a relatively small set of 
"irreducible" integrals for each topology, plus boundary integrals generated when the pf 
factors cancel propagators. Of course the "irreducible" integrals are only irreducible with 
respect to ( [4.16 ), and not with respect to the integration-by-parts and Lorentz identities. 



We compute the "irreducible" integrals once and store them. 

For the planar double box, eq. ( 4.16| ) and momentum conservation, ki-\-k2 + k3 + k4 = 0, 



show that the "irreducible" monomials needed to generate all T^interf. are 

'P^rirti^^ ™) = (2^/ • hr (2p • hr, m + n < 6. (4.17) 

The restriction on the sum of m and n comes from gauge theory — at most six powers 
of the loop momentum can appear in the Feynman diagram numerator algebra. The 
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planar double box also has a symmetry, reflection about the central vertical line, so that 
'Pinterf. ("^) = "Pinterf . "T") • This leaves Only 16 integrals ( 4.17] ) to store. Not all of the 
highest degree integrals actually appear in the amplitude. For the non-planar box, we store 
the integrals of the monomials. 



m + n < 6. 



(4.18) 



We also store those "irreducible" boundary integrals with six (instead of seven) propagators; 
these monomials are generated by three independent factors. For example, for the planar 
double box boundary integral obtained by setting 1/3 — > 0, one requires 



m + n + r < 5. 



(4.19) 



For the helicity amplitude approach, the more general loop momentum polynomial 
Phei. given in eq. (|4.60 requires a bit more work before the above method can be used. 
Consider the product Si ■ p. Because is a four-dimensional vector, this can also be 
written as Si ■ p[4] . We can expand p^^ in terms of a basis of four different four-dimensional 
vectors. Because of momentum conservation, there are only three independent external 
momenta, but we can use the Levi-Civita tensor to construct a fourth one. 



Then 



where 



4 = — 



4 



4 



1 

2su 
1 

2si 
1 

2tu L 
4 



t {2p ■ki)+u {2p ■k2) + s {2p ■ ks) 
t {2p ■ ki) - u{2p ■ k2) + s {2p ■ k^) 
t{2p-ki)+u{2p-k2) - s(2p- fcg) 



(4.20) 
(4.21) 



~ stu ^^^^'^'^^ ~ stu 

Thus we can write 

Ei ■ p = c^Ei ■ ki + c^Ei ■ k2 + ■ ks + c^ei ■ V , 



(4.22) 



(4.23) 



This equation, and the analogous one for Ei ■ q, reduce the problem of handling helicity 
amplitude polynomials ([4.6|) to those of the form 



Vhei. = V{p^,p- q, q^,p- ki, q ■ h; Xj, Ag, A^+g, v-p,v 



(4.24) 



where v ■ p and v ■ q come from the Ct, coefficients. 

The effect of inserting factors oi v ■ p and v ■ q into the integral is very similar to 
inserting factors involving the (— 2e)-dimensional components of the loop momenta, Xp and 
Xq. In either case, shifts of p and q by amounts proportional to the external momenta 
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ki, as required to perform the Gaussian integrals over p and q, have no effect, because 
V ■ ki = \j ■ ki = 0. Thus both types of factors result only in polynomials in the "vacuum 
graph" Schwinger parameters Tq. 

Insertions of A^, and Xp^g may be handled by differentiating the (— 2e)-dimensional 
part of the (Wick rotated) integral, 



d "^^Xp d "^^Xq exp 



X'l rp \'2rp \2 rp 

-'^p ~ -^q -'^q ~ ^p+q ^pq 



with respect to Tp, Tg and Tpg. They lead to parameter insertions such as 



Ap - 


'^q ~l~ '^pq 


^ A 


? 


-^p+q 


Jp+Tq 


^ A 


! 






ApAg - 




T-l2 

-6)^ 
^ A2 ' 


(Ap)^ Aq - 


- e(l - 6) 


^Tq + Tpq 

[ A2 


2 x 2 x 2 

p p+q 


.6(1-6) 


Tp -\- Tq + T^g 


' A2 



T^Tq + 



pqi 



+ (2 



A3 

T T T 

, p-'- q-'- pq 



A3 



Similarly, the polynomials in u • p and w • (7 are easily parametrized: 

StU Tg + Tpg 



{v-pf 
{v-{p + q)f 

{v ■ pf 
{v-pf{v ■ qf 
{v-pf{v-qf 
{v-pf{v ■ qf{v ■ {p + q)f 



StU Tp + Tg 

~8 A ' 

J Stu\^ {Tg+Tpg 



V 8 y 

stu \ ' 

^) 

stu \ ' 

^) 

stu\ ' 



A2 



Tf^ 
A2 



. T„ + Tpg Tpg(Tg + Tj 



A2 

T„ + T, + T 



15- 



pqj 



A3 

T T T 



^ -p ^ ^q^ -^pq _ -^^^p^q^pq 



A2 



A3 



(4.25) 



(4.26) 



(4.27) 



These equations apply to any two-loop integral, independent of the external momenta. 
They also apply in the presence of any additional numerator factor of the form f(p-ki,q-ki), 
since cross contractions are forbidden by the orthogonality of v to the ki. 

Finally, polynomials inv-p, v-q can be related to those in Aj, using the expansion ( 4.21| ). 
For example. 



+ Ap = p[4] • p[4] = s + 1 44 + ^ 44 - T'^*" i4f, 



(4.28) 
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or 



where 



Similarly, 



where 



Vp = —p^ + S C^C2 + t C2C3 + u c^c^ . 



P+1 



stu 



{v-{p + q)f + V, 



pq 1 



(4.29) 
(4.30) 



(4.31) 



pq 



-q + sqcg + iCgCg + nc^Cg , 

<P + qf + Kc? + cDi^ + cl) + t{4 + cl)(c| + 4) + u{4 + cl){4 + 4) . 



(4.32) 



Note that Vp, Vq and "Ppg only contain the types of Lorentz products which already appear 

in Pinterf.- 



Because eqs. ([4.26| ) and ( 4.27 ) are so similar, and taking into account the relations 
between the integrands, we can solve for the additional "new" integrals required for the 
helicity method, in terms of the "old" integrals needed for the interference method. For 
example, for a general function f{p-ki,q-ki),we have 



^l^l f 

\2r„. „^2 



stu 



e {v{p + q)ff 



2e 



1 - 2e 



/ 



e(l-e)/ 8 



3 

4e(l - e 



46(1 - e) 



(1 - 2e)(3 - 2e) 



(1 - 2e)(3 - 2e) 
2 



VpVqf + 



f 



stu 
4 



iv-pnvqrf= — 



1 -2e)(3-2e) 
8e(l-e)(2-6) 



3 - 2e ./ A ' 

e 



+ 



1 - 2e)(3 - 2e)(5 - 2e) 
4e(l-e) fVpf 



VpVqf+, 

'PlVqf + 



2e / A 



(3-2e)(5-2e) y A 



/ 



(4.33) 



A factor of 1 /A indicates that a shift of the dimension of the integral is required: D D+2, 
e ^ e — 1 (e's in prefactors should not be shifted, however). 

In practice, we used eqs. ( 4.29| ) and ( |4.31 ) to eliminate v • p and v ■ q in favor of A^, 
Ag and A^^^ in the loop momentum polynomial. We used equations like ( 4.33| ) to compute 
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the "irreducible" monomials including the Aj, which we then stored. E.g., for the planar 
double box integral, we stored values for 



P, irred 
hel. 



{h,l2, h,m, n) = (Xiy^ (A^)'^ (a2+J'3 (2g . (2p • k^T, 

2li + 2/2 + 2/3 + m + n < 6. 



(4.34) 



Having reduced all the tensor loop integrals in the amplitudes to a linear combination 
of master integrals, the next step is to expand the master integrals in a Laurent series 
in e, beginning at order using results from refs. |53, 54, 47, 46|. Many of these 

master integral expansions are given in terms of Nielsen functions |55|, usually denoted 
by Sn.p{x). However, it is straightforward ||5^ to express the results solely in terms of 
poly logarithms [57|, 



U2{x) 



dt 



Lin-i(t) , 



t) 



(4.35) 
(4.36) 



with n 



2, 3, 4. The analytic properties of the non-planar double box integrals appearing 

; there is no Euclidean region in any of 



in the amplitudes are somewhat intricate |34, 
the three kinematic channels, s, t or u. So we do not attempt to give a crossing-symmetric 
representation, but instead quote all our results in the physical s-channel (s > 0; t, u < 0) 
for the gg 



gg kinematics ( 2.1 



4.2 Checks on results 

We performed a number of consistency checks on the amplitudes to ensure their reliability: 

1. As a check of gauge invariance, we verified that the amplitudes vanish when a gluon 
polarization vector is replaced with a longitudinal one. 

2. The agreement of the explicitly computed infrared divergences with the expected 
form (2.8) provides a stringent check on the amplitudes. Most of the master integrals 
contain divergent as well as finite terms, so the finite remainders are checked indirectly 
in this way. 



Using supersymmetry Ward identities we evaluated the identical-helicity case, 
including fermion loops ||2^, by relating it to the already known identical-helicity 
pure-glue gg gg amplitude [Q. The integration in ref. was done by a com- 
pletely different technique, thus checking the programs and integration methods used 
to obtain the general helicity cases. 

As described in more detail at the end of section ^, we compared our results for 
99 ~^ 99 to those of ref. |12|. The interference of the two-loop gg — > gg helicity 



amplitudes with the tree amplitudes, after summing over all external helicities and 
colors, and accounting for the different schemes used (HV vs. CDR), agrees precisely 



with the calculation of ref. [ 12 1 . 
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4.3 H^'^^ operator 

As mentioned in section ^, the function H^l^{e), which contains only 1/e poles, has not 
been predicted a priori for general processes. However, there is accumulating evidence from 
explicit calculations |ll|, |l^, ^ in the CDR scheme that the color- and helicity-summed 
matrix element (7W(0)|ii'(2)(e)|;\4(0)) 

is a sum of terms for each external colored leg in the 

process, namely 

(A^W|J/(2)(e)|A^W) = 4^^,^K<) + n,/7f))(A^W|A^W) , (4.37) 

where Ug is the number of external gluons, and is the number of external quarks plus 
anti-quarks, with 



9 \ ^AA 19 / \ 79 ins / J AM 07 J ^ ^ ' 



7r2 
'72 


89 \ 

108y 


NN^ - 


Nf 

— + 






--\ 




~ 96 


108 ) 



tt(2) (^^ 11 2 409\ ^2 ( ^'^ 41 \ / 3^ 7r2 3 \ 1 

U""'"^^ + 864j^ +r4^'"^"T^J^V 2^'^^"^J^^ 

+ f^-^l^iV;. (4.39) 
V48 216y iV ^ ^ ' 

(2) (2) 

Note that H), ' and m ' are constants, independent of the kinematic variables. 

We find that the full color and helicity dependence of -ffjjs (e) for gg gg is the sum 
of two terms, 

/ „2 X 2e 

H 

where 



Ti • To , T, • T, 



(4.41) 



with ln((— s)/(— t)) ^ Ins — ln(— t) — ivr in the s-channel, etc. (The overall factor of 
(/i^/(— s))^'^ is a choice of convention, because {jj? / {— s)Y'^ — 1 is of order e. Including it 
cleans up the finite remainder M^gg^gg a bit.) The first term in the sum is proportional 

(2) 

to the identity matrix in both helicity and color spaces. In the HV scheme, Hg is given 



by precisely the same value ( 4.38| ) found in the CDR scheme ||T^. The value in the FDH 
scheme is different; see eq. ( |5.9D . 

The second term in eq. ( [4.40 ) is also independent of the helicity configuration, but 
it is a nontrivial commutator matrix in color space. (The possibility of nontrivial color 
structure in H^'^\e) was pointed out in ref. ||2^.) Indeed, it vanishes when sandwiched 
between tree amplitudes, after performing the color sum, 

oc {M^^'^\(Ti ■ T2T2 ■ - T2 ■ T^Ti ■ T2)\M^^'^) = 0, (4.42) 



using hermiticity of the Tj. Equation ( 4.42 ) ensures that the result (|4.40D is perfectly 
compatible with the previous color-summed results ( [4.37|) . 

Actually, in the course of the qq —>■ qq calculation, the authors of ref. evaluated 
the matrix element 

(m'S^,,Ih'"iM^^,,}, (4.43) 
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relevant for the interference term in identical-quark scattering. Because the color structures 
for the s-channel and i-channel tree amplitudes, M^^^^^iq, and M^^qi^^^i, are different, 
eq. ( [4.42 ) does not apply, and the non- vanishing CDR result is completely consistent with 
our HV result eq. (4.41), including all normalization factors [58]. 

We find that the commutator term is the same in FDH scheme as well. Note that it 
can be rewritten as a "triple product" in color space, because 

Ti • T2 , T2 • Ts] = ifabcT^T^n ■ (4.44) 

This form of the color operator has previously appeared in analysis of the contributions of 
one-loop factors for soft radiation (i.e., S^^^ in eq. ( 3.25| )) at NNLO |26|. This fact, and 

- (2) 

the lack of scheme dependence for H , leads one to suspect that it arises from soft, not 
coUinear, virtual contributions. The 1/e divergence would presumably cancel against the 
contributions discussed in ref. ||2^, in a color-resolved approach to a NNLO computation. 
In a fully color-summed approach, however, such contributions should cancel individually, 
thanks to eq. ( [4.421) . 

The factorization of H in eq. ( 4.41[ ) into a product of kinematic and color factors 
is clearly an accident of having exactly four external colored partons. Color conservation, 
Ti -|- T2 + T3 -|- T4 = in the color-space notation, implies that there are only three 
independent Tj • Tj factors, say Ti • T2, T2 • T3, and Ti • T3. But their sum is a c-number. 



Ti • T2 + T2 • T3 + Ti • Tg 



(4.45) 



hence there is only one independent commutator. (A similar argument holds in the 



fabcT^TjT^ representation [26|.) For three or less external colored partons, all such struc- 



tures vanish; whereas for five or more partons there are multiple independent ones. 



4.4 Finite remainders 

The two-loop finite remainders are defined in eq. (^^) and are color decomposed into 



M 



'A1A2A3A4 ™ Their dependence on the renormalization scale and Nj may 

be extracted as 



M 



(2),[i]fin 
A1A2A3A4 



26o(ln(.//x^)-zvr)Mi;iItA.+ 



M 



(o),H 

AiA2AaA4 



N 



+ 



M 



(2),[i]fin 



/ A1A2A3A4 ^ jY2 ^ A1A2A3A4 ' 
.2\ „-„\ )i^(l).Wfin 



-26o(ln(./A^^)-.vr)Mj^;^^X3A4 + 

Ar2 

_l_ T\T /^[*] _|_ i\T trW J f_ rW 

^ '-^AiA2A3A4 ^ ^ V -"A1A2A3A4 ^ ^AiA2A3A4 ' 



A1A2A3A4 

1,2,3, (4.46) 
7,8,9. (4.47) 



The ^-dependence is a consequence of renormalization group invariance. The tree and 



one-loop functions, -^AiAiA3A4 ^Xi^\2^:!\4^ given in eq. ([232[) and eqs. (|3^ ) - (|3^) , 



r(l).Hfin 



respectively, while 60 and bi are given in eq. (| 
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The coefficient functions A, B, C, D, E, F, G, H, I, which depend only on the Mandel- 
stam variables, obey several relations. Group theory (e.g. U{1) decoupling identities) 
implies that the G functions are not independent of the others, 

^AiA2A3A4 ^(^AiA2A3A4 ~^ ^AiA2A3A4 ~^ ^AiA2A3A4^ ^^1^2\i^A ' (4.48) 



A1A2A3A4 \ A1A2A3A4 A1A2A3A4 A1A2A3A4 / A1A2A3A4 ' V ) 

^AiA2A3A4 ^(^AiA2A3A4 ~^ ^AiA2A3A4 ~^ ^AiA2A3A4^ ^^i^2^z^i ' (4.50) 

and that the sum of the B coefficients vanishes, 

s'^i =-5™ -^['1 f4 5n 

A1A2A3A4 A1A2A3A4 A1A2A3A4 \ ' 

As at one loop, for the ++++, — and 1-+ helicity configurations, Bose 

symmetry under exchange of legs 3 and 4 (t ^ n) implies further relations, 

4''a2A3A4(^'*'^) = 4'i2A3A4(«'^'*)' Xe{A,B,C,D,E,F}, (4.52) 

^ifA2A3A4(^'*'^) = ^i?A2A3A4(^'«'0' Ye{G,HJ]. (4.53) 

In appendix]^ we give the explicit forms for the independent finite remainder functions 
appearing in eqs. ( |4.46D and (4.47). For the two complicated helicity configurations, 1-+ 



and — H h, these functions also involve auxiliary functions, ^sym,[«] ^^^^ j^s-Ym,[i\^ which 

will be presented in appendix |B|. The latter functions also serve as the finite remainders 
for gg ^ gg m N = 1 super- Yang-Mills theory, as discussed in section ^ 

5. Scheme conversion at two loops 

The preceding helicity amplitudes were presented in the HV variant of dimensional regu- 
larization and renormalization. As mentioned in the introduction, the HV scheme contains 
D — 2 = 2 — 2e virtual (unobservable) gluon states, and 2 external (observable) gluon states. 
However, it is possible to alter the number of virtual states. In the FDH scheme |22, 23|, 



one adjusts the number of virtual gluon states to be 2, matching the number of external 
states, and also matching the number of fermionic degrees of freedom in a supersymmetric 
theory. This scheme is quite similar to dimensional reduction (DR) p^ . 

Dimensional reduction is usually thought of as having D < 4, i.e., e > 0, and contains 
D — 2 = 2 — 2e gluon states, plus 2e scalar states, for a total of 2 bosonic states. On the 
other hand, the helicity of a particle is its angular momentum eigenvalue for a rotation in 
the two-dimensional plane normal to its momentum vector. If D is less than four, this plane 
does not exist, making the definition of helicity obscure. The FDH scheme can be regarded 
as an analytic continuation of DR to D > 4, to make it compatible with helicity amplitudes. 
No scalars are required, however. For both the HV and FDH schemes, helicity amplitudes 
with fermions and gluons are computed in the same fashion, with D-dimensional loop 
momenta and four-dimensional gluon polarization vectors (see section [4.1| ). In performing 
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the algebra leading to the loop-momentum polynomial when the trace of the Minkowski 
metric is encountered, one sets 



= Ds=A-2e5R, 



(5.1) 



5r 



1, HV scheme, 
0, FDH scheme. 



(5.2) 



This procedure is gauge invariant because the terms proportional to Dg are related to loops 
containing scalar fields in the adjoint representation |23|. We allow Sr to be arbitrary below, 
although only the HV and FDH cases in eq. ( ^.21 ) seem well motivated. 

The CDR scheme has Dg — 2 = 2 — 2e virtual gluon states, just as in the HV scheme; 
but in addition there are D — 2 = 2 — 2e external gluon states. To convert from the HV 
to the CDR scheme within the helicity method, one could in principle compute additional 
amplitudes where some external states have e-helicities (explicit polarization vectors that 
point into the extra (— 2e)-dimensions) |^9[. Since the CDR result is already available via 
the interference method ||l^, we have not done that computation. Instead we shall check 
the conversion between schemes expected from experience at one loop. 

A given scheme has implications for regularization of both ultraviolet and infrared sin- 
gularities. These implications have been discussed extensively at one loop |^, 33, 6C]. Let 
us first consider the ultraviolet situation. Renormalization by modified minimal subtrac- 
tion, as in eqs. (^) and (|2.6|) , leads to different renormalized coupling constants, related 
by finite shifts. In the class of schemes we are considering, the ultraviolet behavior only 
depends on the number of virtual gluon states. Thus the CDR and HV schemes imply the 
same coupling constant, the standard MS coupling, as{fJ-). The FDH and DR schemes also 
are the same in the ultraviolet (the ability to accommodate helicity, and the sign of e, are 
irrelevant here), and so they both define the DR coupling, a^^{^). 

To shift from either pair of schemes to the other, in the ultraviolet, amounts to using 
the following relations between coupling constants |^, ^ , recently improved to two-loop 
accuracy [p^. 



6 27r 18 V 27r y 



+ 



(5.3) 



6 27r 18 V 27r ' 



(5.4) 



(Recall that the three-loop running coupling enters into any NNLO computation. The 



three- loop beta- function coefficient 62 in DR differs from the value in MS |62|, but it can 
be obtained simply from the coupling shift (|5.3| ) [p^].) For completeness, we give the 
two-loop relation between the MS coupling and that defined by an arbitrary value of 
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6 ^ ^ 27r 



+ 



(5.5) 



6 ('-'^^)^^ + 



C2 



+ 1^(1-^.) 



+ 0([af-(^)]= 



12 



(1 - Sr) 



27r 



+ 



(5.6) 



which reduces to eqs. ( |5.3| ) and (5^) for 5r = 0. Inserting the appropriate couphng relation 
into the perturbative expansion of the amphtude ( p. 4] ) leads to simple, finite ultraviolet 
conversion relations between renormalized amplitudes 

Because the ultraviolet shifts are so simple to implement, in the rest of this paper (with 
one exception to be discussed below) we take the ultraviolet scheme to be the same as the 
infrared scheme. That is, when we report results for ■M.'gg-^gg in the FDH scheme, they 
correspond to coefficients of a perturbative expansion defined as in eq. (|2.4| ), but where 
as(/i) is the FDH/DR coupling, a^^(/i). In the more general 6r scheme, the expansion 
parameter would be af^(/i). 

In this convention, the one-loop relation between gg gg helicity amplitudes is |22, 33| 



\MgLgg)s, = |-M^;,t,,)HV + ^ (1 - Sr) IM^^gLj , (5.7) 
which only involves a finite shift. The corresponding relation at two loops also requires 



+ ^{1-Sr)\M^'^ 
b 



shifts of the divergent terms in the infrared decomposition (2.8). We find that 



Knv-CA[\ + \e]{l-5R), 



Sr), 



fin 



+ 



p2 



TT^ + — 

144 12 



-Cf 



(5.8) 
(5.9) 

Sr) I-m£L,,) , 

(5.10) 



where i^HV is given in eq. ( |2.13| ) and Hg'^'^ (the value in the HV or CDR schemes) is given 
in eq. ( [1.381) . Because I^^\2e,fi; {p}) contains at most 1/e^ poles, the term proportional to 
ex {l — 6ji) in Ksj^ clearly could be shifted into {Hg )sj^ if desired. However, the assignment 
we have chosen makes eq. ( ^.lOD simpler. 
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Also, the interpretation of Kj^ s. as the integral of a splitting function^ ^ leads to 
both the e'^ and e terms proportional to (1 — Sr) in eq. ( [5.8D : The azimuthally-averaged 
g ^ gg splitting kernel is given in the general 5r scheme by 160] 



{P'^{z-e)) = 2CA 



1-z 



+ 



l-z 



+ 1 + 



il-5R)]z{l-z) 



With the identification [60| 



dz iz{l-z))-^ {Pli{z;e)) 



2Ti 



^ + -fa+[Ka--Ti\e, 



(5.11) 



(5.12) 



one sees that the (5R-dependent term oi K = Kg is 



-^{1-6r) Cdz {z{l 
1 - e Jo 



z))^-^ = -CA[\ + le]{l-5R), (5.13) 



in agreement with eq. ( ^.8[ ). 

One can also present results for the two-loop gg gg amplitudes using the Sr scheme 
as the infrared regulator, but switching to the MS coupling constant with the aid of eq. (^). 
For the infrared decomposition (|2.8| ) to hold, assuming that I^^\e) is scheme-independent, 
we find that the quantity -ftTRg, must be set to 



Ksu = Kuv + ^{l-SR) + 0{e) 



(5.14) 



Thus such a definition of K would be scheme dependent too. Also, in contrast to the 
simplicity of eq. (|5.9| ), the scheme-dependent part of iJ"*-^^ will contain logarithms and will 
no longer be proportional to the identity matrix in color space. Hence we refrain from 
presenting such a decomposition explicitly. 

Finally we discuss conversion from the HV scheme results reported in section ^ to the 
CDR scheme used in ref. [^]. In the CDR scheme, one usually computes the interference 
of amplitudes, summed over all external colors and (2 — 2e) polarizations. The generic 
one- loop/tree interference encountered at NLO is 



2Re/<^^°)^2Re [(M^^l-^f) 



(5.15) 



color, hel. 



Inserting the infrared decomposition (|2.7D for A4n^ into eq. ( 5.15| ) gives 



2Re [(-Mfl^^^^lM^ 



(0)\ 



7(l,0)fin 



color, hel. 



where 



(0)\ 



(5.16) 



(5.17) 



color, hel. 



^We thank Henry Wong for clarifying this point. 
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It is well-established from explicit calculations and general arguments |22, 60 1 that 
the finite remainder ( 5.17| ) has the same value in the HV and CDR schemes, in the limit 
e — > 0. Essentially, the treatment of unobserved partons is the same in both schemes, 
so the infrared divergences should take the same form, when expressed in terms of the 
lower-order-in-a^ amplitudes. 

It is natural to expect the same pattern to hold at two loops. The two-loop/tree 
interference is 



color, hel. 



RS 



2Re Y: [(M°^|/(^)|MW> + (M^)|/«t|M°^) 



?(2,0)fin 
~r -'rs : 



(5.18) 
(5.19) 



color, hel. 



where 



r(2,0)fin 

-'rs 



2 Re Yl [{M^n^^''\M 



(o)\ 



color, hel. 



(5.20) 



Note that I^^^ and J^^) 

are the same operators in the HV scheme as in the CDR scheme. 
We have interfered the color-decomposed finite remainders of the two-loop gg gg 
helicity amplitudes in the HV scheme, as given in section ^, with the tree amplitudes given 
in eq. ( 2.22| ), summing over all external helicities and colors with the help of eq. ( 2.24| ). 
This sum gives precisely the same result as the corresponding quantity ( |5.20| ) in the CDR 
scheme, as evaluated in ref. |12|, after accounting for the slightly different definition of H^'^^ 
that we used in eq. (14.401) . We conclude that eq. (|]2^ should be the same in the HV or 
CDR schemes for general two-loop QCD scattering amplitudes. 



6. Two-loop amplitudes in pure = 1 super- Yang-Mills theory 

The quarks of QCD are (massless) fermions transforming in the fundamental representation 
of SU{N). If one replaces the quarks by a gluino, a massless Majorana fermion transforming 
in the adjoint representation, one obtains a supersymmetric theory, pure = 1 super- Yang- 
Mills theory. The amplitudes for this theory, when it is regularized in a supersymmetry 



preserving fashion, obey supersymmetry Ward identities |16], and from experience at one 
loop they are expected to be simpler than the corresponding QCD amplitudes. On the 
other hand, SU{N) group theory generates linear relations between amplitudes of the 
two theories, so one can use the two-loop super- Yang-Mills amplitudes to simplify the 
presentation of the two-loop QCD amplitudes, as we do in appendix 

In this section we discuss the supersymmetry Ward identities and infrared decomposi- 
tion for the two-loop amplitudes in pure = 1 super- Yang-Mills theory. Then we describe 
the finite remainder functions for these amplitudes, deferring the most complicated formu- 
las to appendix p. 

Here we work in the FDH scheme discussed in section ^, in order that the Ward iden- 
tities are valid. One set of identities implies that "maximal helicity violating" amplitudes 
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vanish for any super symmetric theory and any number of loops, 

Mf'^''{gt,gt,9t,---,9i) = 0, (6.1) 
Mf'^''{g^,gt,9t,---,9i) = 0, (6-2) 

where g {g) denotes a gluon (gluino), and the superscripts denote hehcities in the all- 



outgoing convention. In addition to eqs. ( |6.1D and (6.2), all other super- Yang- Mills n-point 
amplitudes containing either zero or one negative-helicity particle vanish trivially, by using 
gluino helicity conservation. We have checked that eq. (|6.l[) is indeed obeyed for the four- 
point amplitude gg — > gg at two loops in the FDH scheme [p^ ]. 

Other identities relate the non-vanishing supersymmetric helicity amplitudes for exter- 
nal gluons alone, to amplitudes where some of the gluons are replaced by gluinos. For the 
four-point amplitudes, all the non-vanishing amplitudes in pure = 1 super- Yang-Mills 
theory can be related to the gg — > gg amplitudes p3[ : 

Mf^{g^rgt.g^,gt) = ^^Mr^ig^,gi,g,,gt), (6.3) 
Mr^iairgtrg^rgt) = ^^Mr^igi,gt,g,,gt). (6.4) 



These relations are crossing symmetric, when a crossing symmetric definition |63] of the 
spinor products is used. Thus, to obtain all the gg gg, gg — > gg, gg — >■ gg, and gg — > gg 
amplitudes from eqs. ( |6.3| ) and ( |6.4| ), it suffices to give the two independent non- vanishing 

helicity amplitudes for gg ^ gg, namely the 1-+ and — i h configurations. 

First we present the infrared decomposition of the pure = 1 super- Yang-Mills am- 
plitudes at two loops. The equations in section || hold with a few modifications for the 
super- Yang-Mills case. We use the perturbative expansion (2.4) but in terms of the FDH 
(or DR) coupling ( |5.3D . The group-theoretic replacements required to convert the quarks 
to gluinos are 

Cf ^ Ca, TrNj ^ Ca/2 . (6.5) 

Some of the previous equations, such as eq. ([4.38|) for Hf'\ are given for gauge group 
SU{N) in terms of N and Nf, with Tr = 1/2, rather than in terms of general Casimir 
operators. In such equations, to recover the Casimir representation, one should first sub- 
stitute Nf 2TRNf and 1/N ^ {Ca - 2Cf), followed by iV ^ Ca. Then one can apply 
the substitutions ( |6.5| ). 

The first two coefficients of the beta function for pure N = 1 super- Yang-Mills theory 

5-- = |Ca, = (6.6) 
and jg is similarly modified, 

Tf" = Ica. (6.7) 

The coefficients K and Hg"^^ are given in the FDH scheme, via eqs. (|5.8| ) and (^.S|), by 

K^S^ = (3 - ^ - ^e) Ca, (6.8) 
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«te=(| + ^-^)cl. (6.9) 



With these replacements, the one- and two-loop infrared decompositions ( p. 7] ) and ( p.q ) 
hold in the super- Yang-Mills case. We also color decompose the amplitudes and strip off 
the helicity phases exactly as in eq. ( ^.15 ) for the QCD case. The dependence of the one- 



and two-loop finite remainders on ^ and N are then extracted as 

<SSr = ^^A^fk ' ^ = 7>8,9, (6.10) 



M 



(2),SYM,[i]fin 
Ai A2A3A4 



^f-)2 (ln(.//x^) - iTrr + (In(VM') " i^) 



M 



26-- (ln(.^^) - zvr) M«^;™,f + iVM 



(o),H 

Ai A2A3A4 

SYM,[jl 



A1A2A3A4 '^AiA2A3A4 ' 

i = 1,2,3, (6.11) 



M 



(2),SYM,[i]fin 
Ai A2A3A4 



(ln(./;x^) - ^n) M^^.^^^'^ + A^G^-'S,^ , 



i = 7,8,9. (6.12) 

The one-loop supersymmetric remainder functions are given in terms of the QCD ones, 

i = l,2,3, (6.13) 
i = 7,8,9. (6.14) 



a 



A1A2A3A4 



SYM,[lJ 

^Ai A2 A3A4 



A1A2A3A4 A1A2A3A4 g A1A2A3A4 ' 

2/'qSym,[1] sym,[2] , ^SYm,[3] \ 

\^ a1a2a3a4 a1a2a3a4 a1a2a3a4 j 



The M^'^'^'W correction term in eq. ( |6.13 ) is a consequence of using two different schemes — 
FDH for super- Yang-Mills theory vs. HV for QCD. The two-loop analogues of eq. (p^ ) 
are eqs. (|A3^ ), ( |06D , ([OtI ), ( [AlgD , ( [AISOD , (|A]5|), ( [EsBD , and (|A]5|) in appendix 
These equations also have correction terms due to the different schemes used, as well as 
feed-down from the subtracted singular terms, which depend on the fermion representation. 
The correction terms are more complicated than at one loop, involving Af^^^'W^" as well as 
j\^{o),W^ g^ill they contain no special functions, only logarithms. 

The two- loop supersymmetric remainder functions A^™, B^^^, G^^^^ obey the same 
types of identities as the corresponding QCD functions. The group theory relations are. 
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SYM,[7] 
A1A2A3A4 



+ A 



,m \ 
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Ai A2A3A4 



2(^SYM,[1] ^SYM,[2] 

A1A2A3A4 A1A2A3A4 
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(6.15) 
(6.16) 
(6.17) 

(6.18) 



The Bose symmetry relations, which hold only for the 



-^^SYM,[2]/ , N 



X 



SYM,[1] 



s,u,t) , 



-++ helicity configuration, are 
X^{A,B}. (6.19) 
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In appendix we give the explicit forms for the independent = 1 supersymmetric 
finite remainder functions appearing in eqs. ( |6.11 ) and ( |6.12| ). This completes the de- 
scription of the two-loop four-point amplitudes for N = 1 super- Yang-Mills theory, and 
simultaneously of the auxiliary functions required for the QCD amplitudes presented in 
section ^. 



7. Conclusions 

In this paper we have presented the two-loop amplitudes for gluon-gluon scattering in 
QCD, and for all of the 2^2 scattering processes in pure N = 1 super- Yang-Mills theory, 
including the full dependence on external colors and helicities. We found that there is an 



r(2) 



additional 1/e pole term, H in eq. ( 4.41 ), which has nontrivial color dependence, and 
which vanishes after interfering it with the tree amplitude and summing over colors. We 
investigated the dependence of the amplitudes on the flavor of dimensional regularization 
employed. The QCD results, when summed over all external colors and helicities and con- 
verted to the CDR scheme, are in complete agreement with the previous results of Glover, 
Oleari, and Tejeda-Yeomans |12|. We also expressed the one-loop-squared contribution to 
the NNLO gg — > gg cross section in terms of one- loop finite remainders. Again the ap- 
propriate interference, converted to CDR scheme, is in complete agreement with previous 
resuhs g. 

Much numerical work still remains in order to implement the two-loop amplitudes of 



this paper, or those of refs. |11, 12], in a numerical program for NNLO jet production at 
hadron colliders. When that is accomplished, however, the intrinsic precision on the QCD 
predictions should reach the few percent level, providing a stringent test of the Standard 
Model at short distances. 
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A. Finite remainder functions for QCD 



In this appendix, we present the explicit forms for the independent finite remainder func- 
tions for gg — > gg in QCD, which appear in eqs. ( 4.46| ) and ( 4.47 ). For the helicity 
configuration, the functions are 
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B. Finite remainder functions for pure = 1 super- Yang-Mills theory 

In this appendix, we present the independent = 1 supersymmetric finite remainder 
functions appearing in eqs. ( 6.11 ) and ( 6.12| ). For these functions are 
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